Synchronization on complex networks with different sorts of communities 
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In this paper, inspired by the idea that many real networks are composed by different sorts of 
communities, we investigate the synchronization property of oscillators on such networks. We iden- 
tify the communities by the intrinsic frequencies probability density g{u}) of Kuramoto oscillators. 
That is to say, communities in different sorts are functional different. For a network containing two 
sorts of communities, when the community strength is strong, only the oscillators in the same com- 
munity synchronize. With the weakening of the community strength, an interesting phenomenon. 
Community Grouping, appears: although the global synchronization is not achieved, oscillators in 
the same sort of communities will synchronize. Global synchronization will appear with the further 
reducing of the community strength, and the oscillators will rotate around the average frequency. 

PACS numbers: 89.75.Hc, 05.45.Xt 



I. INTRODUCTION 

Many social, physical and biological systems are often 
described by networks. The collective dynamical behav- 
iors of those networks are simultaneously determined by 
the individual systems, the interactions between individ- 
ual systems, and the network structure. Very recently, 
it is found that the structure of many large-scale real 
networks are neither regular nor random, but of some 
common and important characters, such as short average 
distance, large clustering coefficient and the power-law 
degree distribution [ij]. Networks of such structure are 
called complex networks. Besides the characters men- 
tioned above, many real-world networks have the so- 
called community structure 0, 0] ■ community in a net- 
work is a set of nodes, where the edges inside are much 
denser than those connecting the nodes belong to this 
set and the rest of the network. Examples of community 
networks are numerous in biological, technical and social 
systems. It is common in society that people are divided 
into different groups by their sex, age, interests and so 
forth, and each group can be considered as a community. 

One of the most significant aims of the studies on com- 
plex networks is to understand the effects of the network 
structure on the dynamical processes, among which the 
synchronization of coupled oscillators is the simplest but 
one of the most important ones. Synchronization phe- 
nomenon has been observed for hundreds of years and 
in a variety of field, including natural, physical, chemical 
and biological systems Because of the limitation of 
knowledge on network structure, the studies of synchro- 
nization are restricted to either on the regular lattices or 
on the random networks for a long time period. Recently, 
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with the pioneer works on small- world [5| and scale-free 
network models @ , much attention has been paid to the 
studies of synchronization on those complex networks. 
Soon, scientist found delightedly that oscillators on com- 
plex networks are much easier to synchronize than on 
regular lattice with the same size and density of edges 
0,11, , which further stimulates their enthusiasm. 

Up to now, the relationship between network struc- 
ture and synchronizability [uf, [H [H, [li, [11, [H, ^ is 
well understood and many ways to enhance the network 
synchronizability [H, [H, HO, |2l|, [H [H [13, [H, p, [13, 
Hlj [29I [30I have been proposed (see also the review ar- 
ticle [si'l and the references therein). Very recently, the 
synchronization properties of networks with community 
structure are investigated, it is found that the connect- 
ing pattern between communities have great effects on 
the synchronization process |33i] , and the community 
structure will hinder the global synchronization of oscil- 
lators [H, [3^ : the stronger the community structure the 
worse the global synchronizability. All the works men- 
tioned above concentrate on the networks with only one 
sort of communities, i.e., all the communities are similar 
in both function and topology. However, in real world, 
many networks have communities with different sorts, an 
example in point is a friendship network of children in a 
U.S. school Hi]. As shown in Fig. 3.4 of Ref. [s^, 
there are four communities: white middle school stu- 
dents, white high school students, black middle school 
students, and black high school students. Clearly, the 
four communities can be divided into two sorts by race: 
white and black, and there are two communities in each 
sort. Or, those communities can be divided by age, as 
middle school and high school. 

In this paper, with the help of Kuramoto model 
[37I [ssl . [sol . [40I [llj . we investigate the synchronization 
properties of complex networks with different sorts of 
communities. In our study, oscillators on different sorts 
of communities are identified by the natural frequencies 
probability density g{Lij), namely, the intrinsic frequen- 




FIG. 1: (Color online) The relationship between order pa- 
rameter M and the coupling strength a for the network with 
two sorts of communities, and the number of communities in 
those two sorts are equal. Red circle, green triangle, blue up 
triangle and light blue diamond curves represent the order 
parameters of Ai, Bi, A and B, respectively. And the black 
square curve represents the order parameter of the whole net- 
work. Each point is obtained from 100 independent runs. 

cies of nodes on different sorts of communities are taken 
differently. 

This paper is organized as follows. In section II, the 
Kuramoto model and the order parameter are intro- 
duced. In section III, the network model we investigate 
will be described in detail. And also, we will give the sim- 
ulation results of synchronization properties of Kuramoto 
oscillators on complex networks with different sorts of 
communities. The conclusion remarks are drawn in sec- 
tion IV. 



FIG. 2: (Color online) The averaged sine of phases of os- 
cillators vs. time for community networks with C = 0.3 at 
different coupling strengths. There are two sorts with the 
same number of communities in the network. 
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II. KURAMOTO MODEL AND ORDER 
PARAMETER 

In this paper, we use the coupled phase oscillators, 
Kuramoto model [13, [H, [H, SO, |4l| , to analyze the col- 
lective synchronization on complex networks. A modified 
Kuramoto model can be described by the coupled differ- 
ential equation (ssj : 



(1) 



where 0^, uji are the phase and the intrinsic frequency 
of node i, Ai is i's neighbor set, and a is the overall 
coupling strength. The intrinsic frequency Ui is chosen 
from a probability density g{ui). 

Usually, the intrinsic frequencies of different oscillators 
are assigned differently, if there are no couplings, all the 
oscillators will rotate independently, thus the phases of 
the oscillators are distributed almost uniformly in the 



FIG. 3: (Color online) The relationship between order pa- 
rameter M and the coupling strength a for the network with 
two sorts of communities. Sort A and sort B contain 3 and 
7 communities, respectively. The inset shows the change of 
order parameter in a larger scale. 



interval [0, 27r]. With the increasing of the coupling 
strength, the oscillators will adhere to each other to some 
extent. At last, when the coupling strength reaches some 
critical point CTc, collective synchronization of all oscilla- 
tors emerges spontaneously, although the intrinsic fre- 
quencies of different oscillators are different. 

To measure the synchronization phenomena, an order 
parameter M is introduced: 



M = 



1 



N' 



(2) 
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FIG. 4: The relationship between actual frequencies fli and 
the intrinsic frequencies uji for the network, which is consisted 
of two sorts with equal number of communities. 

where (• • •) and [• • •] denote the average over time and 
over different configurations, respectively. N' is the num- 
ber of nodes that are taken into account. Clearly, AI is 
of order 1 / ^/N' if the oscillators are completely uncou- 
pled (cr = 0), and will approach 1 if they are all in the 
same phase. In our simulation, the dynamical equations 
integrate using the Runge-Kutta method with step size 
0.01. The order parameters are averaged over 2000 time 
steps, excluding the former 2000 time steps to allow for 
the relaxation to a steady state. Note that, we consider 
not only the order parameter of all the oscillators in the 
network, but also investigate the synchronization of par- 
tial oscillators, such as the nodes belong to an individual 
community, or the nodes belong to a sort of communi- 
ties. For different cases, the sum in Eq. (2) goes over the 
oscillators that are taken into account, and N' is taken 
accordingly. 



III. NETWORK MODEL AND SIMULATION 
RESULTS 

Our numerical simulations are based on the commu- 
nity network model in Refs. [s^ll^]- The model starts 
from n community cores, each core contains toq fully con- 
nected nodes. Initially, there are no connections among 
different community cores. At each time step, there are 
n new nodes being added, each node will attach m edges 
to the existing nodes within the same community core, 
and simultaneously m' edges to the existing nodes out- 
side this community core. The former are internal edges, 
and the latter are external edges. Similar to the evolu- 
tionary mechanism of Barabasi- Albert networks we 
assume the probability connecting to an existing node i 
is proportional to z's degree A;,. Each community core 
will finally become a single community of size Nc, and 



FIG. 5: The relationship between actual frequencies Hi and 
the intrinsic frequencies uii for the network, which is consisted 
of two sorts with different numbers of communities: sort A 
contains 3 communities, while sort B contains 7. 



the network size is = nNc- By using the rate equa- 
tion (similar to the analytical approaches used in Refs. 

Hj]), one can easily obtain the degree distribution of 
the whole network, p{k) oc k~^. We simply take the pro- 
portion of the external edges, C = , , to measure the 
community strength. Clearly, a smaller C corresponds to 
sparser external edges thus a stronger community struc- 
ture. 

In our simulation, a network of 10 communities is taken 
into account, and each community contains 200 nodes, 
thus there are totally 2000 nodes in the network. The 
average degree is set as fc = 6. The communities are 
divided into two sorts A and B, and each sort contains 
5 communities, which are labeled as Ai, • • • , A5 and Bi, 
B5, respectively. In fact, all the ten communities 
have similar structure, but the intrinsic frequencies of 
nodes in A are assigned randomly in [0, 1], while those in 
B are assigned in [1,2]. 

To get the synchronization property of this kind of 
networks, we investigate how the order parameters in 
three scales, a community, a sort of communities and 
the whole network, change with the coupling strength. 
Figure 1 shows the simulation results. From top to bot- 
tom, C increases from 0.01, 0.1 to 0.3, corresponding to 
weaker and weaker community structure. In the case 
C = 0.01, with the increasing of coupling strength, the 
order parameter of individual community soon reaches 1 
(the red circle and the green triangle curves), indicating 
the synchronized state within one community. However, 
the phases of oscillators in different communities are dif- 
ferent, even for those belong to the same sort (the blue 
up triangle and light blue diamond lines), and the order 
parameter for the whole network is much less than 1 (the 
black square line). The order parameter for each com- 
munity is near to 1, for community sorts A and B, the 
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order parameter are almost equal and much lower than 
those for individual community, and the order parame- 
ter of the whole network is the lowest one. In the case 
C = 0.1, the oscillators in the same community will also 
synchronize soon. With the further increasing of cou- 
pling strength, an interesting phenomenon appears: os- 
cillators belong to the same sort are synchronized before 
the appearance of global collective consensus. Although 
the communities are equal in structure, the intrinsic fre- 
quencies divide the communities into two groups. We 
call this phenomenon Community Grouping. In the case 
C = 0.3, when the couphng strength is not too large, 
the Community Grouping can also be observed, while for 
sufficiently large cr, this phenomenon disappears and the 
global synchronization is achieved. In figure 2, the av- 
eraged sine of phases of oscillators are plotted with time 
at steady states for community networks with C — 0.3. 
When a — 0.5, the coupling is too weak to form any syn- 
chronization. When a = 3.5, Community Grouping phe- 
nomenon appears. When the coupling strength is strong 
enough (cr = 9.0), all the communities rotate coherently. 
From figures 1 and 2 we know that only in some proper 
regions of (C, cr), the Community Grouping phenomenon 
could appear. We also investigate the network with more 
than two community sorts, and find that if the numbers 
of communities contained by one sort are equal, the net- 
work exhibits similar synchronization properties to the 
case of two sorts. 

As to networks with two sorts containing different 
numbers of communities, some interesting phenomena 
appears. Figure 3 shows the changing of the order pa- 
rameters with the coupling strength, where sort A and 
sort B contain 3 and 7 communities, respectively. For dif- 
ferent community strengths C = 0.01, 0.1 and 0.3, order 
parameters for individual community all reach 1 soon, 
showing that oscillators within one community synchro- 
nize quickly. For the networks with very strong commu- 
nity structure (C = 0.01), even the communities belong 
to the same sort could not get synchronized each other. 
Meanwhile, for larger C, the order parameters for sort A 
and B are not equal: with the coupling strength increas- 
ing from 0, both Ma and Mb rise soon, but the former 
is faster. While, with the further increasing of the cou- 
pling strength, the rising of the two order parameter slow 
down, and the rising of Mb becomes a little faster, sooner 
or later Mb will surpass Ma- This crossing behavior of 
Ma and Mb can be observed in the panel with C = 0.1, 
and larger C will make the crossing point smaller. Af- 
ter the crossing point. Ma will increase, decrease and in- 
crease again, and finally reaches 1 (see the inset of Figure 
3) . The abnormal change of Ma can be simply explained 
as follows. When a is small, with fewer oscillators, A will 
show better coherent. For larger cr, interactions between 
oscillators are strong, because of containing more com- 
munity, B has stronger power of influence, thus oscillators 
in A arc disturbed by B. This effect could enhance the 
global synchronization, but simultaneously reduce Ma- 

When there are couplings between oscillators, the os- 
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FIG. 6: (Color online) The states of the communities. The 
circle and the square present two sort of community and the 
communities with the same color indicate they are synchro- 
nized, (a), (b) and (c) represent the cases of C = 0.01, 0.1 
and 0.3, as shown in the above numerical simulations. 



cillators will not rotate at their own frequencies uj but at 
actual frequencies il, and the distribution of depends 
on the coupling strength and the structure of the net- 
works. Figure 4 and 5 shows the relationships between 
actual frequencies and intrinsic frequencies of each os- 
cillator at different coupling strength given C ~ 0.3. In 
figure 4, there are two sorts with equal number of commu- 
nities. It can be seen that for a weaker coupling strength, 
the actual frequencies ili are almost equal to the intrin- 
sic frequencies w^, increasing the coupling strength a will 
make the actual frequencies of oscillators on the same sort 
more consensus, till the oscillators rotate at a common 
frequency near their own average frequencies, indicated 
by the two parallel lines in the third panel. Further in- 
creasing the coupling strength will make the two parallel 
lines nearer and nearer, and finally, all the oscillators ro- 
tate at about the average frequency of the whole system, 
indicating the achievement of the global synchronization. 
Figure 5 reports the results of network consisted of two 
sorts with different numbers of communities. The actual 
frequencies will change alike, and the final synchronized 
frequency is about the average frequency 1.2. 



IV. CONCLUSION AND DISCUSSION 

In conclusion, we investigated the synchronization 
properties in complex networks with different sorts of 
communities and found that when the community struc- 
ture is strong, only the oscillators in the same commu- 
nity synchronize, with the weakening of the community 
strength, oscillators in the same sort of communities will 
synchronize independently, only when the community 
structure is not evident that all the oscillators in the net- 
work can synchronize, this is clearly shown in figure 6. 
When the numbers of communities in different sorts are 
not equal, the sort with fewer communities shows bet- 
ter consensus at weaker coupling and when the couplings 
between oscillators are strong, the sort with more com- 
munities have larger order parameter. When all the os- 
cillators are synchronized, they will rotate at the average 
frequency. The huge-size real-life networks are generally 
consisted of many communities, and those communities 
may be functional different. In this sense, to study the 
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detailed dynamical properties of network with function- 
ally different communities is of interesting. The current 
work provides a start point on this issue, which is helpful 
for the in-depth understanding of synchronization pro- 
cess on networks. 
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